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ABSTRACT 

We estimate the baryon mass density of the Universe due to the stars in galax- 
ies and the hot gas in clusters and groups of galaxies. The galaxy contribution is 
computed by using the Efstathiou, Ellis & Peterson luminosity function, together 
with van der Marel and Persic & Salucci's mass-to-light versus luminosity relation- 
ships. We find f2|* ars ~ 0.002. For clusters and groups we use the Edge et al. X -ray 
luminosity function, and Edge & Stewart and Kriss, Cioffi & Canizares' (gas mass 
)-luminosity relations. We find Q? as ~ 0.001. The total amount of visible baryons is 
then fib — 0.003, i.e. less than 10 per cent of the lower limit predicted by standard pri- 
mordial nucleosynthesis, implying that the great majority of baryons in the Universe 
are unseen. 

Key words: Key words: galaxies: fundamental parameters - intergalactic medium 
-cosmology: observations - dark matter. 



1 INTRODUCTION 

The estimation of the baryon mass density of the Universe 
essentially involves, for each class of objects having a visible 
baryon content, an integration over luminosity of the prod- 
uct of the luminosity function (LF), 4>(L), the luminosity, 
L, and the mass-to-light ratio for the baryon component, 
Mb/L, according to the expression 



Pb 



(i) 



scales as h% Q (see equation 1), when estimated by dynamical 
arguments the baryon density parameter Sib = pb/pc with 
p c — 5 ftfo x 10~ 30 g cm 3 , the critical density is independent 
of h 5 o- 

In practice, the detailed information required by equa- 
tion (1) has not previously been available. Standard esti- 
mates of pb have assumed a typical value for the visible mass- 
to-light ratio of galaxies, usually inferred from the observed 
dynamics and supposed to be representative for galaxies of 
all luminosities and Hubble types. Then equation (1) reduces 
to 



where T represents E/S0 galaxies, spiral galaxies, clusters 
and superclusters, and (j>(L) and Mb/L refer to the relevant 
class of objects. In the case of galaxies, equation (1) deserves 
some comments. First, let us note that the baryon mass-to- 
light ratio corresponds to the mass-to- light ratio of stars, and 
not to the dynamical one which usually includes dark mat- 
ter (DM). In fact, the presence of DM in the internal parts 
of spiral galaxies is very evident (Persic & Salucci 1991), so 
we cannot ignore its contribution to the mass of a galaxy. In 
addition, the importance of DM is a function of luminosity: 
low-luminosity galaxies are affected more strongly than high 
luminosity ones (Persic & Salucci 1988, 1990). Therefore, it 
is essential that the mass-to-light ratio of the visible baryon 
content be the quantity that enters equation (1). Further- 
more, we emphasize that different populations are present 
in galaxies of different Hubble types. These have different 
LFs, different mean Mj, / L ratios, and different scaling prop- 
erties of (Mb/L) with luminosity. Finally, note that, since pb 



Pb — C 



(2) 



where £ is the galaxy luminosity density obtained by 
integrating the galaxy LF over luminosity and < Mb / L > is 
an assumed mass-to-light ratio. In addition, the hot gas in 
clusters and groups of galaxies has often been neglected in 
the computation of fib 1 

For purposes of illustration, in Table 1 we summarize 
previous calculations of Qb according to equation (2). Note 
that in most cases the dynamical rather than the stellar 
mass to- light ratio is used. The range in the adopted values 
found in the literature does not reflect observational uncer- 
tainties but, rather, real differences in stellar populations, 
in proportions of DM and in the reference radius where the 



1 Throughout the paper, /150 
Ho is the Hubble constant. 



Ho/(50kms 1 Mpc 1 ), where 
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Table 1. Notes/") Adopted luminosity density, in 10 8 h 50 L Q Mpc- 3 (from W Peebles 1971 l 2 l Shapiro 1971; ' 3 1 Kirshner, Oemler & 
Schechter 1979; M Schechter 1976; ' 5 1 unspecified; 16J ^ Efstathiou et al. 1988; M unspecified)/ 6 ' Adopted mass-to-light ratio, in solar 
units. ( c ) Type of adopted mass-to- light ratio, classified according to its procedure of determination. 



Author £(°) < M/L ><» TypeM Q b 



Peebles 1971 


1.5W 


10 h 50 


dynamical 


0.02 


Gott et al. 1974 


O.5I 2 ] 


(2.5 - 7) hso 


dynamical 


> 0.001 


Olive et al. 1981 


2 PI 


(4 - 10) h 60 


dynamical 


0.003-0.007 


Boerner 1988 


0.5W 


(7 - 10) /iso 


4.8 dynamical 


0.011-0.016 


Hogan 1990 


1.51 5 ! 


2 


unspecified 


0.008/igo 1 


White 1990 


l[6] 


1 5 


stellar 


0.003/1" 1 


Kolb & Turner 1990 


1.2 PI 


< 5 h 50 


dynamical 


< 0.01 



dynamical mass-to-light ratio is considered. Detailed infor- 
mation has recently been published on the properties of stel- 
lar mass-to-light versus luminosity relations for E/S0 and S 
galaxies (van der Marel 1991; Persic & Salucci 1991), as 
well as on the relation between hot gas content and X-ray 
luminosity for clusters of galaxies (Edge & Stewart 1991a, 
b). Also, LFs have become available for E/S0 and S galax- 
ies separately (Efstathiou, Ellis & Peterson 1988) as well 
as for clusters (Edge et al. 1990). It is therefore possible to 
obtain dynamical measures of the visible mass associated 
with these structures. The aim of this paper is to use these 
advances in estimating the value of fit,. The value we find, 
fit, ~ 0.003 is lower than some previous estimates. Com- 
pared to current standard nucleosynthesis predictions, our 
estimate aggravates the problem of the 'missing baryons'. 



2.1 Ellipticals 

From Efstathiou et al.'s (1988) field-galaxy LF we get <f>, — 
8.5 x 10- 4 hl Mpc- 3 , a = 0.48 and L„ = 3.3 x 10 10 /i,r 2 L Q ; 
we take x m in — 0.04 and x max = 8. There is strong evidence 
that ellipticals do not have DM inside their optical radii, 
so their dynamical masses at the optical radius and their 
stellar masses coincide (see Djorgovski & Davis 1987). From 
van der Marel (1991) we have A — 4}i 5 oMq/Lq and rj — 
0.35 (all values refer to the B band). Inserting these values 
into equation (5), the integral takes a value of 0.97, and 
we then get p b = 7.4 x 10~ 33 /i5 g cm~ z . Thus the baryon 
contribution of elliptical galaxies to the mean density is 

n (B+so) = 1 5 x 1Q -3 ( E + g0 galaxies ) ( 6 ) 



2 THE MEAN BARYON DENSITY FROM 
GALAXIES AND CLUSTERS 

In all cases, the LF will have the usual Schechter (1976) 
form: 




where (j>„ is a normalization constant, a is the slope of the LF 
at low luminosities, and L* is the luminosity corresponding 
to the 'knee' of the LF; and the baryon mass-to-light versus 
luminosity relation will be a power law: 

with A the mass-to-light ratio at the characteristic luminos- 
ity L*. Inserting equations (3) and (4) into equation (2) we 
get 

p b = <f>*L*A I x L - a+v e- x dx (5) 

where x = L/L* and x min and x ma x represent the ob- 
served minimum and maximum luminosities for a given class 
of objects. The values of the parameters </>*, a, L*, x min , 
x m ax, A and T are observationally known and will be cho- 
sen accordingly for each class of objects. We shall treat sep- 
arately elliptical (and SO) galaxies, spiral galaxies, groups 
and clusters of galaxies. 



2.2 Spirals 

From Efstathiou et al. (1988) we get 0* = 1.1 x 
10~ 3 /i!o MpcT 3 , a = 1.24 and L, = 4.2 x lO lo /i^ o 2L ; we 
take x m i n = 0.01 and x m ax = 8. Dark and visible matter 
are already well mixed in the optical regions of spirals (e.g. 
Persic & Salucci 1991). From Persic & Salucci (1990) dynam- 
ical disk/halo mass decomposition of rotation curves, we get 
A — I.2/15OM0/L0 and rj = 0.35 (all values refer to the B 
band). Inserting these values into equation (5), the integral 
takes a value of 0.94, so we get p b — 3.5 x 10~ 33 /i§o 9 cm~ 3 . 
corresponding to 

Q ( b S> ~ 0.7 x 10~ 3 (S galaxies). (7) 

2.3 Rich clusters 

The hot, X-ray emitting, diffuse intracluster gas is a distinc- 
tive prominent baryon component of clusters (e.g. Sarazin 
1986). From Edge et al. (1990) we get = 6.9 x 
10~ 8 hlo Mpc~' A , a = 1.65 and L, = 8.1 x 10 44 h M f erg/sec; 
we take x m in = 0.012 and x ma x = 2.5; relative to the 
gas mass within 0.5/i^ Mpc radius, from Edge & Stewart 
(1991a) we take A = 50 h^ B g erg' 1 sec and rj = -0.62 
(All the X-ray data are relative to the (2 - 10) keV band.) 
Inserting the above values into equation (5), the integral 
takes a value of 7.6, corresponding to (a contribution to) 
Of, = 1.4 x 10~ 4 . This estimate refers to the gas mass con- 
tained within 0.5h^Q Mpc radius (i.e., within - 2 optical 
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core radii, see Bahcall 1977). In, at least some cases, how- 
ever, there is evidence (Jones & Forman 1984) for significant 
amounts of gas at larger distances, ~ 3h^ ( l Mpc. This dis- 
tance is a good measure of the typical virialization radius of 
rich clusters [e.g., from a recent study of a large rich-cluster 
sample we have R vir = (3.0 ± 0.2)h^ Mpc (Biviano et al. 
1992)], and corresponds to the Abell (1958) radius, ie. the 
radius roughly encompassing most of the cluster's member 
galaxies. The hot-gas surface-brightness profile can be de- 
scribed in terms of the so-called /3-model, 

S(r) = S [l + (r/a) 2 ]-^ +1/2 

with a (the core radius) and f3 as free parameters, which 
implies a de-projected spatial gas density distribution of the 
form 

Pga s(r)=p [l + (r/a) 2 }- 3 ^ 2 

(e.g., Henriksen & Mushatzky 1985; David et al. 1990). Fits 
to Einstein IPC data far a survey of clusters give for /3 an 
average value of < f3 >= 0.66 ± 0.10 (see Jones & Forman 
1984; David et al 1990). Therefore, on average, beyond a 
few core radii the hot- gas density falls off Pgas T , 
and the mass in gas rises linearly with radius. We therefore 
estimate that the hot-gas contribution of rich clusters from 
within the Abell radius is 

^(ci.gas) ~ 8.6 x W~ 4 h^ B (intracluster gas) (8) 

In addition to the diffuse gas, further baryons in clusters are 
contributed by the stellar component of the member galax- 
ies. However, recent direct estimates (see Edge & Stewart 
1991b) indicate that, in the cores of clusters, the stellar mass 
is up to a factor of 3 smaller than the gas mass. Further out 
there are indications that the galaxy distribution drops more 
rapidly than the gas (e.g., Eyles et al. 1991), so that within 
a few core radii the stellar mass is significantly lower than 
the gas mass (roughly 10 per cent for a stellar mass-to-light 
ratio of 4; see David et al. 1990), and at the Abell radius the 
stellar component is negligible compared to the gas mass. 

2.4 Poor clusters and groups 

A population of poor clusters and groups of galaxies must 
certainly contribute to the baryon content. The information 
on the LF and the gas content at these mass scales is quite 
scanty. We extrapolate the Edge et al. (1990) cluster X- 
ray LF through the range of luminosities typical far these 
structures, 41 < logL(2 - lOfceV) < 43 (see Kriss, Cioffi & 
Canizares 1983; Bahcall, Harris & Rood 1984; see Bahcall 
1979 for a discussion of the continuity of cluster and group 
LFs). Thus in equation (5) we take the parameters of the 
Edge et al. (1990) cluster LF (see above), with x m i n — 1.2 x 
10~ 4 and x m ax = 0.012. Relative to the hot-gas content 
within 0.5/i^q 1 Mpc radius, the Einstein data of Kriss et al. 
(1983) imply A = 65h^ h5 and n = -0.5 (in the 0.5-4.5 keV 
band). The integral in equation (5) takes a value of 13, so 
we derive pb = 1.6 x 10~ 33 g cm -3 h!~l{f corresponding to a 
contribution to Q,b — 3.2 x 10~ 4 hgg ' 5 Allowing far any hot 
gas extending out to a typical group virialization radius ( 
~ lhV^Mpc, see Pisani 1990), by the same argument used 
for the rich clusters we get 

^(gr.gas) _ 6 4 x 10 -* h -i-8 (intragroup gas). (9) 



3 DISCUSSION 

From equations (6)- (9), we conclude that the combined 
baryon contribution of galaxies and clusters to the mean 
density is: 

Q. b ~ 2.2 x 10~ 3 + 1.5 x lO" 3 /^ 1 - 5 (10) 

As 1 < /150 < 2 with a currently favored value of = 1.5 
(e.g., Pierce & Tully 1988), the cluster/group density contri- 
bution(see equations 8 and 9) is probably Q, gaa — 0.001.. The 
estimated contribution of the stellar populations of galax- 
ies (see equations 6 and 7) is fi* ~ 0.002. This value is in 
very good agreement with the cosmological mass density of 
the damped Lya system, which are un likely to be proto- 
galactic discs and would give a present-day contribution of 
Sl b (Lya) ~ 0.002(see Wolfe 1988). Note that the explicit 
inclusion of the general trend of decreasing stellar mass-to- 
light ratio with decreasing luminosity and the accounting for 
different morphological classes have the effect of reducing 
the estimated baryon contribution of galaxies. In fact, the 
often-assumed values of Mb/ Lb = 8/150 for ellipticals and 
Mb/ Lb = 3^50 for spirals actually refer only to the bright- 
est objects (i.e. Lb > 6L*). Assuming such values as typical 
would overestimate the baryon contribution from stars by a 
factor of - 2. As an example (of it), let us take our mass-to- 
light ratios at L* as the typical values (i.e. 4 far ellipticals 
and 1.2 far spirals) to be used in equation (2). Let us also 
use separate LFs to compute the luminosity densities due 
to ellipticals and spirals (i.e., 0.25 x 10 8 Lb q hV^Mpc^ 3 and 
0.54 x 1O s Lb h^Q Mpc~ 3 respectively, from Efstathiou et al. 
1988). Then equation (2) would be written as 

p b = [C < {Mb I Lb) >]e/so + {£< (Mb/Lg) >] s 

= ([0.25 x 10 8 x 4] + [0.54 x 10 s x 1.2])h^ M Q /Mpc~ 3 

Thus we would find (from the above equation) contributions 
to fib, galaxies of 1.4 x 10~ 3 and 0.9 x 10 -3 for the two classes 
separately. This example shows that the discrepancy be- 
tween some previous estimates of fi» (claimed to agree with 
nucleosynthesis limits) and our own estimate arises mainly 
because these calculations used dynamical and unrealistic 
mass-to- light ratios and, to a lesser extent, because they 
did not use separate LFs for different morphologies. Further 
baryons, in the farm of hot (T > W 7 K) diffuse gas, could be 
supplied by superclusters. However, the observational limit 
to any diffuse X-ray emission from candidate supercluster 
cores (see Persic et al. 1990), alongside the spatial frequency 
of superclusters (ie., ~ 10 8 ft 3 Mpc~ 3 ), indicates a negligible 
contribution: 

Q SC,gas _ 10 -B 

Another reservoir of baryons could be cold diffuse HI 
gas in the local intergalactic medium (IGM). However, the 
absence of the Gunn-Peterson (1965) trough on local scales 
places a severe upper bound to any such contribution: 
Qhi < 5 x 10~ 7 (See Davidsen et al 1977) 

We therefore expect that our derived value of Q b in 
equation (10) gives a close representation of the actual vis- 
ible baryon content of the local Universe. [Fabian's (1991) 
recent estimate of the baryon density in the Shapley su- 
percluster, fij, ~ 0.18 h^Q 5 is not necessarily in contra- 
diction with our estimate. In fact, Fabian's high value re- 
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lates to a most exceptional region of ~ 35 h§ Mpc size 
and with Sp/p ~ O.Sfi^ 1 , while our low value refers to the 
whole nearby Universe out to a radius of ~ 300h^ Mpc 
and with Sp/p ~ 0] Our value, fl b ~ 0.003, is substan- 
tially lower than the prediction of standard cosmic nucle- 
osynthesis, 0.04 < flbhl with a probable value of n& ~ 0.06 
(e.g., Kolb & Turner 1990; Peebles et al. 1991). Comparing 
the observed and the predicted baryon abundances, we con- 
clude that the stars and gas of galaxies and clusters/groups 
account for only < 10 per cent of the primordially synthe- 
sized baryons (see also Hogan 1990; Kolb & Turner 1990). 
Where are the 90 per cent of missing baryons? We have only 
taken an inventory of the visible baryons associated with vis- 
ible structures. Additional baryons, unaccounted for by the 
present census, can either be clumped in some dark form, 
for instance forming dark haloes around galaxies, or be dis- 
tributed in a diffuse ionized background. [Note that the dif- 
fuse DM associated with clusters, exceeding the gas mass 
by at most a factor of - 3 (e.g., Eyles et al. 1991), could not 
solve the discrepancy with standard nucleosynthesis even if 
it were completely baryonic] Let us consider each of the two 
possibilities in turn. 

(i) Baryons in haloes. The dark haloes of spiral galaxies 
may well extend out to 10-20 times the size of the optical 
discs. Thus it is not difficult to conceive that, by integrating 
the dynamical mass-to-light ratio of galaxies (computed at 
such extended radii) over the LP, the nucleosynthesis value 
of fij, ~ 0.06 might be easily reached. An attractive way 
of hiding the missing baryons, therefore, is to assume that 
they constitute the DM in galaxy haloes. This possibility 
may find support from the evidence that cooling flows may 
be producing baryonic DM in the form of low-mass stars 
or brown dwarfs (e.g., Fabian, Nulsen & Canizares 1984). 
Based on cooling flow analogies, Thomas & Fabian (1990) 
have argued that baryonic DM only forms on fairly large 
mass-scales where gas cooling is quasi-static, while Ashman 
(1990) has suggested that baryonic DM forms on galactic 
and subgalactic scales, following rapid gas cooling (Ashman 
& Carr 1991). 

(ii) ionized IGM. The second possibility is to suppose 
that galaxy formation is extremely inefficient, so that only 10 
per cent of gas in the Universe is now in collapsed structures 
such as galaxies. This could arise if gas was never incorpo- 
rated into galaxies, or if gas was expelled from protogalaxics 
by supernova explosions or galactic winds (Bookbinder et al. 
1980). In this scenario, most baryons in the Universe now 
constitute a smooth ionized IGM, in agreement with the lack 
of any Gunn-Peterson trough in the spectra of quasars both 
at high redshifts (see Sargent & Steidel 1990) and locally 
(Davidsen et al. 1977). There is no shortage of proposed 
sources of ionization. Among conventional ad hoc sources 
there are Population III stars (Carr 1990), dwarf and star- 
bursts galaxies (Silk, Wyse & Shields 1987; Songaila, Cowie 
& Lilly 1990) and obscured quasars (Miralda-Escude' & Os- 
triker 1990). Alternatively, a non-conventional self consis- 
tent explanation involves the radiative decay of cosmologi- 
cal dark matter particles (Rephaeli & Szalay 1981; Sciama 
1982). We conclude by stressing that the result presented in 
this paper, in connection with the dynamical estimates of the 
total mass density of the Universe on scales ~ Wh^Q Mpc, 
Qo > 0.2 (sec Davis et al. 1980; Davis & Peebles 1983) sup- 



ports the conventional view on the need for non-baryonic 
extra dark matter. 
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